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Abstract 



' This paper concerns multiphase models of tumor growth in interaction with a surround- 

O I- ing tissue, taking into account also the interplay with diffusible nutrients feeding the cells. 

Models specialize in nonlinear systems of possibly degenerate parabolic equations, which in- 
clude phenomenological terms related to specific cell functions. The paper discusses general 
modeling guidelines for such terms, as well as for initial and boundary conditions, aiming at 
both biological consistency and mathematical robustness of the resulting problems. Partic- 
ularly, it addresses some qualitative properties such as a priori nonnegativity, boundedness, 
and uniqueness of the solutions. Existence of the solutions is studied in the one-dimensional 
time-independent case. 
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1 Mixture-theory equations for tumor growth 

q [ 1.1 Mixture-theory-based models 

The interest toward mathematical modeling of tumor growth rose considerably in the last decades, 
to such an extent that it has now become one of the most studied topics in mathematical biology. 
Early mathematical models [15, 16, 17] considered tumors as ensembles of only one type of cells. 
Growth was described under the main assumption of constant cell density, by relating the volume 
variation of the tumor mass to birth and death of cells triggered by nutrient supply. In most 
cases simple in vitro geometries were considered, such as spheroids, and qualitative analyses of 
the resulting free boundary problems were detailed [12, 13, 20, 21, 25, 26, 27]. 

However, the biological literature pointed out soon that tumors should be regarded more prop- 
erly as ensembles of different interacting components, e.g., normal and abnormal cells, intercellular 
fluid, extracellular matrix. This aspect is taken into account by modeling tumors as multiphase 
materials by methods of mixture theory, see for instance [14]. In mixture theory one introduces 
a few volume ratios (f> a , where the index a labels the components of the mixture, expressing the 
percent amount of the constituents. Each volume ratio is supposed to satisfy < (j> a < 1, with a 
possible further condition <j) a = 1, called saturation constraint, if one assumes that no voids 
are left within the mixture. Mass balance equations are written for the constituents under the 
assumption of same density: 

^ + v-(0 Q v Q ) = r a , (i) 

ot 

where v Q and T a are the velocity and the source/sink term of the constituent a, respectively. 
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A first class of multiphase models is obtained from Eq. 1 via suitable closure relations relating 
the velocities to the volume ratios of the constituents [1]. A very common assumption is that 
all cell populations share the same velocity, while non-cellular components have their own. Then 
geometrical considerations, for instance some symmetries, may determine kincmatically the ve- 
locities of some constituents, as it happens in [6] where the cylindrical symmetry of tumor cords 
developing along a blood vessel is used to deduce the velocity of the cells. The model presented 
in that paper combines ideas coming from mixture theory with free boundary issues. Particularly, 
it includes dynamical constraints on the nutrient distribution across the tumor mass for modeling 
the formation of an outer necrotic shell of dead cells at the periphery of the tumor cord. 

A second class of multiphase models is obtained by joining to Eq. 1 some stress balance 
equations for the constituents of the mixture, in which inertial effects are neglected: 

-V- (</> a T a ) + a Vp = m a . (2) 

Here T a , m a are the excess stress tensor and the resultant of the external actions on the con- 
stituent a, respectively, whereas p is the intercellular fluid pressure. Equations 2 are used to derive 
the velocities v Q from mechanical reasonings on the internal and external stress sustained by the 
constituents. For instance, assuming that the external actions m Q can be expressed as viscous 
frictions, thus involving only the relative velocity of pairs of constituents, one gets (possibly gen- 
eralized) Darcy's laws that can be plugged into Eq. 1. Considering in particular a sub- mixture of 
two cell populations, namely tumor cells labeled with a = T and healthy host cells labeled with 
a = H, which share the same mechanical properties, and assuming that the extracellular matrix 
only acts as a rigid non- remodeling scaffold providing them with a support for their movement, 
the following equations are obtained [32]: 

= r a , (3) 

where <j> := </>t + 4>h is the overall volume ratio of the cellular matter, E = is the intercellular 
stress (such that Tt = Th — — £(</>)I), and finally K Qm is the motility tensor of the cell population 
a within the extracellular matrix. 

The source/sink terms r Q model proliferation or death of cells, taking into account both natural 
processes linked to the vital cell cycle and the availability of nutrient. Hence they depend on the 
cell volume ratios <pT, <Ph and on the nutrient concentration c, which entails T a — T a {4>T, 4>h, c). 
This introduces the new variable c, for which an evolution equation, usually of reaction-diffusion 
type, is supplied: 

J-V-(DVc)= ]T Q a (4> a ,c), (4) 

a=T, H 

where D is the diffusivity tensor and Q a models the absorption of nutrient by the cells of the 
population a. The tensor D may be taken independent of c, and also of the cell volume ratios <j) a , 
because nutrient molecules are not regarded as a part of the mixture, i.e., they are assumed to 
diffuse through the mixture without occupying space. Equation 4 features then a linear diffusion. 
In addition, the functions Q a are sometimes linear in c, hence Eq. 4 turns out to be, in most 
cases, a linear model for the nutrient concentration. However, the simultaneous dependence of the 
functions r Q , Q a on both <j> a , c makes Eqs. 3, 4 ultimately coupled. 

Technical simplifications in Eqs. 3, 4 may involve the assumption of homogeneous isotropic 
motility of the cells in the extracellular matrix, as well as homogeneous isotropic diffusion of the 
nutrient through the mixture. These imply K Qm = n am l in Eq. 3 and D = DI in Eq. 4, for 
positive constants n ami D. In addition, invoking the hypothesis of same mechanical properties 
for tumor and host cells, one may set urm = ^Hm ='■ K m > 0, which, as pointed out in [32], is 
a good approximation at least in the initial stages of the development of a tumor, when contact 
inhibition among cells is more influential than differences in their motility. 



dt 
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Figure 1: The domain ft in case of cell segregation. At each time, tumor cells are contained in 
f2r(i) and healthy host cells in The two cell populations are separated by the interface 

S(t), which is a {d — l)-dimensional manifold in f2. 



1.2 Cell segregation 

An interesting case is when tumor and host cells remain segregated, i.e., the spatial domain 
il C R d of the problem can be split at each time in two (open) sub-domains Q,T(t), Slij(i) such 
that 0, T (t) U fiff(t) = ft, and H 0, H (t) = (Fig. 1). In practice, each sub-domain contains a 

mixture of extracellular fluid, extracellular matrix, and just one type of cells obeying the following 
balance equation: 

^ - KmV • [&,V(^ a £(^a))] = T ^ a , c), (5) 

each volume ratio <j) a being defined only in the corresponding domain Q a . Equation 5 is derived 
from Eq. 3 noticing that <fi = <p a in Q a owing to segregation. The two mixtures interact at the 
interface S(t) := dQ.T{t) n dVLnit) separating the sub-domains, hence each <j) a solves in principle a 
free boundary problem because S(t) is not fixed. However, it is possible to supplement Eq. 5 with 
proper conditions on S(t) so as to reformulate it globally in £1 Proceeding in a formal fashion, 
we integrate Eq. 5 on f2 Q (i) up to a certain final time T max > 0, then we apply Gauss' Theorem 
to the divergence term at the left-hand side and sum over a to discover: 



E 

a=T, H 



^ Tmax ^max ^ 

/ / ~tf dX dt ~ Km J J 4>oN {4>o^{(t> a )) ■ ndadt 
\ o n a (t) o dn a (t)\s(t) J 



+ n m j j [(t>TV(<f>TX(<t> T ))-<l>HV(<f>HX((l>H))}- ndadt 



S(t) 



E 



r Q (0 Q , c) dxdt, 



a = T , H o n a (t) 

where da is the (d — l)-dimensional Hausdorff measure in R d and n is the outward normal unit 
vector to the boundary on which integration is performed. In particular, along S{t) it denotes 
the outward normal unit vector to so that the analogous vector for Qr(t) is — n (but, of 

course, the opposite convention may also be adopted). 

Next we reintroduce the function (f> : [0, Tmax] xSJ^I: 

\<f)H{t,X) ltX€ilH{t), 
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and use U a Cl a (t) — 0, U a d£l a (t) \ S(t) — dil to obtain 

T max T max 

y J^dxdt-n m J j ' 4>V '{4>T.{(j>)) ■ nda dt 
on o 90 

T max T m 



where we have defined 



+ K m J J (<£))]■ ndcrdt = y j T{t, x, (f), c)dxdt, (6) 

o s(t) o n 

r(i, s, 4>,c):= T <*{4> a , c)l na{t) (x). (7) 



In formulas 6, 7, [•] denotes jump across S(t) whereas lfj a (t) i s the indicator function of the set 
0, a (t). If we reapply Gauss' Theorem to the second term at the left-hand side, we can regard Eq. 6 
as the integral version of the differential equation 

^ _ KmV . [0 V (0E(0))] = r(t, x, c) (8) 

provided n m \<jN (<j>Y>(cf)))\ ■ n = on S(t) at each time. With this condition, Eq. 8 is equivalent to 
Eq. 5 on either sub-domain 0, a (t), being at the same time posed globally in fl We will come back 
later (cf. Sect. 3.2) to the significance of such interface condition from the modeling viewpoint. 

By comparing Eq. 8 with the standard mass balance equation of continuum mechanics we infer 
that the velocity v of the cellular matter is 

v = -/t ro V(^E(0)). (9) 

Furthermore, by defining 

$'(«) := *(«£(«))' (10) 
we notice that Eq. 8 can be rewritten in the form of nonlinear diffusion: 

^-K m A*(<f>) = r(t,x,<t>,c). (n) 

Conversely, Eq. 4 is naturally defined on the whole il: segregation is for cells, not for nutrient. 
Nevertheless, coherently with the segregation assumption, we redefine the right-hand side as 

Q(t, x, <j), c) := Qa((t>, c)ln a (t)(z) 

a=T,H 

and write 

-£-DAc = Q(t,x,<t>,c). (12) 



1.3 Aims and scope 

This paper is concerned with mathematical models of tumor growth of the kind outlined above, 
with a twofold goal. On the one hand, to discuss biologically consistent modeling lines for the 
phcnomcnological terms of the equations (namely, the functions E, T, and Q), as well as suitable 
boundary, interface, and initial conditions. On the other hand, to obtain qualitative results, such 
as a priori nonnegativity, boundedness, and uniqueness of the solution, along with continuous 
dependence estimates, which support modeling with mathematical rigor. For this reason, the 
paper is ideally divided in two parts. 

The first part, encompassing Sects. 2, 3, is especially devoted to modeling. Particularly, Sect. 2 
surveys the most popular models proposed in the literature for E, T, and Q. Inspired by them, it 
fixes some modeling assumptions which will be used throughout the subsequent sections. Section 3 
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discusses boundary, interface, and initial conditions needed to formulate mathematical problems, 
with special emphasis on the use of the former for simulating the surrounding environment e.g., a 
nearby vasculature. 

The second part, encompassing Sects. 4-6, is targeted at analytical issues. Specifically, Sect. 4 
is a preliminary technical one, introducing the main notations and recalling the essential theo- 
retical background. Subsequently, Sects. 5, 6 approach the time-dependent and time-independent 
problems, respectively, establishing a priori estimates on their solutions. Existence of solutions is 
also explicitly addressed in the one-dimensional stationary case. 

Finally, Sect. 7 sketches some research perspectives on recent multiphase models of tumor 
growth incorporating explicitly the attachment /detachment of cells to/from the extracellular ma- 
trix. 

The paper is equipped with two Appendices, which contribute to make it as self-contained as 
possible. Appendix A concerns the handling of the nonlinearities in the equations of the models. 
Appendix B further extends the theory to other kinds of boundary conditions, partly different 
from those discussed in Sect. 3, relevant for applications. 

2 Constitutive assumptions 
2.1 The cell stress function 

The function E expresses the internal stress to each cell population. As already mentioned, the 
Cauchy excess stress tensors are given by Tt = Tjj = — E(</>)I, hence E acts as an intercellular 
pressure depending on the local cell packing. 

For theoretic purposes, in the sequel it will be more customary to deal with the so-called 
constitutive function defined by Eq. 10, rather than with E itself, although it is obviously 
possible to switch at any time to either function via the above-mentioned relationship. 

Diffusion problems are well-known to be ill-posed if the diffusion coefficient is negative, there- 
fore a very basic requirement in our case is that $' be nonnegative. A more complete characteri- 
zation is provided by the following assumption: 

(HI) $ : M — > R is smooth, strictly increasing, and normalized in such a way that $(0) =0. 

Strict monotonicity of the constitutive function is classically required in the theory of nonlinear 
parabolic equations [37]. Notice that if $ is strictly increasing then $' cannot vanish but possibly 
at the origin (indeed $'(0) = is forced by Eq. 10 if (sE(s))' is not infinite in s = 0), therefore 
$'(s) > for all s ^ and $ is invertible. We anticipate that we will use invertibility in Sect. 6 
for the existence theory of the solutions to the stationary problem. 

Many models in the literature assume that E grows steeply as cells get highly packed. For 
example, in [35] a particular instance of the following function is proposed (Fig. 2, left): 



where a, b are positive constants with b >• a, n > 1 is integer, cfi* e (0, 1) is the close-packing 
cell volume ratio, and (•)+ denotes the positive part of its argument. In practice, such a E is 
a physiological pressure for normally packed cells, which rapidly increases as soon as the tissue 
becomes overly dense. Function 13 is smooth for n > 1 and piecewise smooth for n = 1, and the 
corresponding constitutive function $: 



E(s) = as + b (s-0») + 



(13) 



$(s) = -as 3 + b {s -(/)*) 



+ 



n 




s{s-^) + [(s-0*) + f 



n+1 (n + l)(n + 2) 



fulfills hypothesis (HI). 

Other authors use [30, 31] (Fig. 



2, center) 



E(s) 
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Figure 2: From left to right, the intercellular stress functions 13, 14, 15. 

A, t > 0, which grows again linearly for small s (physiological pressure) but then reproduces a 
release of the stress after the maximum r/(2v / A) attained for s = l/\/X. This should model a 
saturation effect due to that at high densities all cells are not able to push. It is interesting to 
note that such a behavior is opposite to the one assumed by Eq. 13, nevertheless from Eq. 10 it 
can be easily computed 

r , , t I arctan (y/Xs) s \ 

m = -x{-^=x TTa?J' 

which complies in turn with hypothesis (HI). 

In the previous two examples it results E(s) > for all s > 0, but this is not strictly necessary 
for hypothesis (HI) to be satisfied: even if E(s) < for some s > one might get $'(s) > for all 
s 7^ 0. As recalled in [1, 14], negative values of the stress model adhesive intercellular forces, which 
compete with the repulsive ones because cells, unless too packed, like to stick together to form 
multicellular aggregates. This behavior is reproduced, for instance, by the function E proposed in 
[34] (Fig. 2, right): 

s 



-log 

s 



if n = 1 

S 0* 

-^_>' ifn>l, 
. n — 1 s 

where now 0» denotes the stress-free volume ratio corresponding to unstressed tissue (i.e., £(</>*) = 
0). Notice that E(s) < for \s\ < <fi*, with E(s) — > — oo for \s\ — > 0. The resulting constitutive 
function: 

$(«) = Isl""^ 

is strictly increasing for all n > 1 and turns Eq. 11 into the porous medium equation with nonlinear 
forcing term. Function 15 somehow summarizes qualitatively the trends of the functions 13, 14 at 
large volume ratios, indeed for 1 < n < 2 it is bounded from above and, if n < 2, tends to zero 
when s — > +oo (saturation effect), whereas for n > 2 it grows unboundedly and more and more 
steeply as n increases. 

In [22] the authors introduce the idea, next borrowed by a few other papers [1, 14, 19], that 
the cell stress blows up when </> approaches a maximum allowed volume ratio (/> max <= (0, 1], which 
corresponds to an asymptote of the function E for s — <f) max (Fig 3, left): 

S(«)=p(0ma»-0*) | y (16) 

|s|(0max - S) 

where p > is a constant coefficient and 0* denotes again the stress-free volume ratio. Notice 
that, E being infinite at </> ma x, is also infinite, which violates hypothesis (HI). However, 
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^ ^* ^max ^ ^* ^1 ^2 ^max 

Figure 3: From left to right, the intercellular stress function 16 and that used in [14]. 

strict monotonicity of $ on (-co, </> max ) i s preserved by function 16, while in general it fails on 
[0, 0max) with the function proposed in [1, 14] (Fig. 3, right). The latter is such that S(s) < 
for < s < <p2 with a local minimum at s = <px e (0*, (f> 2 ) in order to take into account cell 
adhesiveness at low volume ratios, then E(s) > for (f> 2 < s < </> max with S(s) — > +oo when 
s — > max to reproduce cell repulsion. In addition, they set E(s) = for < s < <f>* to model 
that far apart cells ignore each other. We refer the reader to [14] for the analytical expression of 
such a S. In [19] the adhesion region is instead eliminated by taking = 4> 2 and setting £ = 
before 0*. In practice, the resulting function is the positive part of function 16, which restores the 
monotonicity (however not strict) of $ before max . 

It can be argued that condition E(s) — > +oo for s — > </> max is mainly intended to enforce the 
bound (j) < max on the solution to Eq. 8, trusting to the physical intuition that a strong, infinite 
in the limit, intercellular pressure prevents cells from packing too much. Nevertheless, we will 
prove that this is not necessary to achieve the proper upper bound on <j). 

Before concluding the discussion on the constitutive function, we introduce a condition that 
will play a fundamental role in the forthcoming theory (cf. also [24, 29]): 

Definition 1. We say that a function / : K — > R is &-Lipschitz continuous in an interval ICR, 
with constant Lip $ (/) > 0, if 

|/(* 2 ) - /(si)| 2 < Lip*(/) ($(s 2 ) - $(*i)) {s 2 - si), Vai, s 2 g /. 
Notice that the right-hand side is nonnegative because $ is increasing. 



2.2 The growth term 

The function T expresses proliferation or death of cells in connection with the availability of 
nutrient. The basic principles inspiring the modeling of T are usually that few cells proliferate less 
than many cells, that proliferation stops when cells fill all the available space, and that cells need 
a minimum amount of nutrient to survive. 
For instance, in [11] they use 

m C ) = (i- )_^_-^_^_, 

where the first term at the right-hand side is the cell growth due to mitosis (fostered by the 
availability of nutrient), the second term is the cell death (enhanced by the lack of nutrient), and 
So, . . . , S4 > are parameters. Notice that cell growth is zero whenever = (no cells) or 
<t> = 1 (cells occupy the whole space, recall the saturation constraint) and that nutrient chemistry 
reminds of the Michaelis-Mentcn kinetic. 
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A simpler form of T, including an explicit nutrient threshold triggering the switch between cell 
proliferation and death, is that proposed in [34] : 

T(ct>, c)=7#l-Mc-c) 

where 7, c* > are parameters. In this case, for <j> <G [0, 1] it results r>0ifc>c*,r<0 
if c < c», hence T acts as a source or a sink, respectively, according to the values taken by c. 
Unfortunately, such a T has the drawback of vanishing for <f> = 1 regardless of c, which implies 
that cells no longer die once they have reached their maximum concentration, even if the nutrient 
falls below the survival threshold. This difficulty may be overcome by the following correction: 

T(<j>, c) = 7^(1 - 4>){c - c*)+ - 7 20(c - c.)-, (17) 

where (-) + , (•)" denote positive and negative parts of their arguments and 71, 72 > are param- 
eters. 

In [3, 4] proliferation and death of cells are linked to energy reasonings, specifically ATP con- 
sumption through nutrient oxidation along the cell cycle, and the following form of T is proposed: 

m c) = ^<f>(f(4>)9(c) - 6)+ §^(/(<^( c ) - §)-, (18) 

where k, 9, Q® M , Ti/ 2 > are the reaction rate of oxidation, the total rate of ATP consumption, 
the average cost of the full cell cycle, and the half-life of dying cells, respectively. The function 
/ (of 4> alone) is introduced to inhibit proliferation (namely, ATP production) in overly dense 
tissues, whereas the function g (of c alone) defines the effectiveness of the oxidative process in 
terms of nutrient supply. In more detail, / vanishes when <fr attains the maximum threshold </> ma x 
and g increases with c. Examples of such functions are f{4>) — <f) maK — <j>, g(c) — c. 

In [19, 32] it is suggested that cell duplication and death may occur on a stress-induced basis, 
depending on the level of compression felt from the surrounding tissue: 

F a (4>, c) = la cf>H{Y,* a - S(0)) (j- - lj - S^HQty) - £**), (19) 

H being (possibly a mollification of) the Heaviside function and S* < £** two stress thresholds. 
When the actual stress acting on the cells is above the first threshold proliferation is inhibited. 
If the stress further grows above the second threshold then cell apoptosis is triggered. Since the 
sensitivity to the stress affects the way in which a cell runs through its vital cycle, which is what 
mainly breaks down when mutations change a normal cell into an abnormal one, the thresholds 
S* , £** are expected to be different for tumor and host cells, in particular S^*- 1 > S^*' [19]. 
Inspired by the examples above, we see that a convenient structure of the growth terms is 

r a (0,c)= J2 -f a m)9 v a {c)-s<t>, 

v—p, d 

which allows one to account for possible differences in the mechanisms of proliferation (y = p) and 
death {v = d) of tumor and host cells, as depicted for instance by Eqs. 17, 19. In more detail, 
the coefficients 7^ are specific proliferation and death rates for tumor and host cells. The terms 
/a( ( / ) )5a( c ) r efer to joint stress-nutrient induced proliferation and apoptosis, under the assumption 
that the cell stress is directly determined by the volume ratio <fi. Finally, the term —5<ft accounts 
for natural cell apoptosis without influence from the distribution of nutrient, the coefficient 5 being 
the same for both cell populations. 

Some technical requirements on the previous terms are now stated, taking into account a 
generic maximum volume ratio </> max € (0, 1] allowed for cell packing 1 : 

0max < 1 then 1 — (/>max is the constant volume ratio of the rigid non-remodeling extracellular matrix, and 
<Pe '■= 0max — <j> the volume ratio of the extracellular fluid (not explicitly modeled) filling the interstices within the 
mixture to enforce the saturation constraint. 
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(H2) 7 g, 6 > 0, 7^ < 

(H3) bounded, $-Lipschitz continuous, and nonnegative in [0, </> max ] 

(H3.1) fP > in (-oo, 0), fg < in (0 max , +oo), /*(0) = /£(</> max ) = 
(H3.2) H < in (-oo, 0), f a > in (0 max , +oo), /*(0) = 
(H3.3) nondccreasing in [0, </> max ] 

(H4) locally bounded, Lipschitz continuous, and nonnegative in R. 

Although not explicitly required, we incidentally notice that further assumptions may be sug- 
gested by physical considerations. For instance, nutrient-induced proliferation may be non- 
decreasing and nutrient-induced death g% non-increasing. 

2.3 The absorption term 

The function Q describes the uptake of nutrient by cells. A common and simple prototype of this 
term is (see e.g., [2, 14, 34]) 

Q{4>, c) = -A0c, 

where A > is a parameter. This form translates the basic principle that the absorption of 
nutrient depends simultaneously on the number of cells present in the domain and on the quantity 
of nutrient available to them, including that few cells uptake, on the whole, few nutrient even if 
the latter is abundantly supplied, and, conversely, that few nutrient can poorly feed a large cell 
population. More general absorption terms are introduced in the series of papers [6, 7, 8, 10]: 

Q{4>, c) = ~H>(c), 

where ip is then suggested to be of Michaelis-Menten type (cf. also [11, 30]). This allows the 
authors to generalize Q to the case of multiple species of nutrients, including in <p the dependence 
on the various concentrations [9]. 

In [4] the absorption term is directly linked to the chemical mechanisms internal to the cells 
responsible for proliferation and death, and the following form of Q is proposed: 

Q{<f>, c) = -\<j>f(cj))g(c), 

where A > is the oxygen uptake rate and the functions /, g are the same as in Eq. 18. 

In view of these examples, and of possible differences in the consumption of nutrient by normal 
and abnormal cells, we refer to the following structure of the absorption term: 

Q a (<t>, c) = -\ a h a (4>)q a (c). 

Here, \ a is the specific absorption rate of the cell population a while the functions h a account for 
cell-dependent uptake dynamics, which may differ for cancer and host cells because of different 
internal chemistry. Finally, q a is the chemical consumption rate of nutrient, which instead is much 
likely to be the same for tumor and host cells as it depends essentially on the chemical properties 
of the environment and of the nutrient itself rather than on cell genetics. 
Some technical assumptions on these terms are in order, namely: 

(H5) X a > 

(H6) h a $-Lipschitz continuous and nonnegative in K 

(H7) q a , locally bounded and nonnegative in [0, +oo) 

(H7.1) qr a (0) = 0, q a < in (-oo, 0) 
(H7.2) q a nondecreasing in [0, +oo). 

Again, additional assumptions, not strictly needed for theoretic issues, may be welcome for 
physical consistency. For instance, one may require the cell-dependent absorption rate to vanish 
when no cells are present, i.e., h a (0) — 0. We anticipate that we will actually use this assumption 
when addressing the existence of stationary solutions. 
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3 Boundary, interface, and initial conditions 



3.1 Boundary conditions 

The parabolic character of Eqs. 11, 12 calls for conditions on the whole boundary dfl x (0, T max ] 
in order to properly define the corresponding mathematical problems. The most common bound- 
ary conditions in tumor growth problems refer to characteristic values of some quantities at the 
periphery of the cell tissue (Dirichlet conditions) or to their fluxes across the external shell of the 
portion of tissue under consideration (Neumann or Robin conditions). Often boundary conditions 
account for interactions of the cell aggregate with the outer environment (not explicitly modeled). 

Let us consider, at first, the chemicals nourishing the cells. Nutrient is supplied to the cells 
by the external environment e.g., by a nearby vasculature, whence, dissolved in the extracellular 
fluid, it diffuses through the cell tissue. The presence of the vasculature can be modeled as a 
boundary condition for Eq. 12: one identifies a portion of the boundary, say dbfl the subscript 
b standing for "blood", with a blood vessel, whence the nutrient carried by blood flows into the 
tissue. One has therefore a condition on the normal flux — Z?Vc • n, n being the outward normal 
unit vector to db^l- This may be either a Neumann condition, if the flux is given, or a Robin 
condition, if the flux is expressed in terms of c itself, for instance as n(c — ct) where rj, Cb > are 
parameters. Such a condition implies that the flow of nutrient across the vessel wall depends on 
the quantity of nutrient already present in the tissue, with respect to a characteristic threshold Cb- 
In more detail, if the concentration c equals the characteristic concentration Cb then there is no 
flux, while if c < Cb then the nutrient flows from the blood to the tissue. Conversely, if c > q, then 
the nutrient flows from the tissue to the blood, namely blood carries away the nutrient in excess. 
However, this situation is not expected to happen because blood is the only source of nutrient in 
this model. The parameter Cb can be identified with the average physiological concentration of 
nutrient in the blood, whereas rj is a characteristic speed of filtration through the vessel wall. We 
point out that Robin's is the biologically most appropriate condition to be imposed at the vessel 
wall, see [28, 36], and in this respect Neumann condition should be regarded as a zeroth-order 
approximation. For this reason, we ultimately set 

-DVc • n = n(c - c b ) on d b £l x (0, T max ]. 

In the remaining portion of the boundary, say df£l the subscript / standing for "far" (in the 
sense that this boundary is far from the vasculature where tumor growth mainly takes place) , one 
might prescribe the concentration c& for conveying the idea that the quantity of nutrient is there 
the average physiological one in healthy tissues: 

c = c b on dfCt x (0, T max ], 

which is a Dirichlet boundary condition. 

We suppose db$l U <9/f2 = <9f2, Int c^Q fl Int = for consistency (cf. Fig. 1) but we do not 
exclude that either dfQ or <9f,f2 is empty, if for instance the vasculature fully surrounds the cellular 
tissue (dfQ = 0, dQ = 9(,0) or if one is concerned with avascular tumors = 0, dtt = dffl). 

Concerning the cells, since we are considering in situ tumor growth, it is reasonable to assume 
that they do not penetrate the vasculature, which corresponds to no flux across the boundary dbfl: 

K m V$(^) • n = on d b fl x (0, T ma J. 

Recalling Eq. 9, this condition can be rewritten as <pv ■ n = 0, which says that the normal 
component of the velocity of the cells vanishes at the blood vessel. Cells neither cross the vessel 
wall nor detach from it, but they can slide along the vessel because no restriction is imposed on the 
tangential component of their velocity. Instead, at the far boundary <9/fi the tissue is essentially 
relaxed, for the main dynamics is expected to occur near the vasculature. Therefore, considering 
that in the present setting the stress state is determined by the cell volume ratio, the natural 
condition for cf) on dffl is some physiological value 

<i> = e (0, max ) on dffl x (0, 
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Figure 4: Left: blood vessels and far boundaries are among the edges of 17. Right: blood vessels 
are internal holes to 17, whereas the far boundary coincides with the outer edges of 17. 

possibly coinciding with the stress-free value if the function £ admits one. 

Boundary conditions outlined here are mainly indicative, and may be detailed more precisely 
for specific domains 17. In particular, the boundaries <9t,17, <9/17 need not be connected, especially 
when several blood vessels are present (in which case <9;,17 will be presumably the union of several 
connected components). As an example, one may consider the applications illustrated in [32], 
which deal with two-dimensional tumor growth in a rectangular domain 17. In some cases, blood 
vessels coincide with one or more (not necessarily adjacent) edges of 17, the remaining ones forming 
instead d/17 (Fig. 4, left). In other cases, all of the outer edges of 17 define <9/17, whereas blood 
vessels are circular holes within the rectangle. The union of their circumferences is then <9;,17, 
which plays the role of an internal boundary to 17 (Fig. 4, right). 

3.2 Interface conditions 

The interface S(t) separating the two sub-domains 17r(i), 17_jf(£) is a material one for the cells, 
meaning that the latter cannot detach from it on either side. This entails the continuity of their 
normal velocity across S(t), i.e., recalling Eq. 9, 

[vnJ = - Km [V(0£(<£))]-n = G, Vt G (0, T max \. (20) 

In addition, classical theory of continuum mechanics requires the continuity of the stress of the 
interfacing materials: 

[</>£(0)n]=O, V<G(0, T max ]. 

For a continuous stress function S, this is fulfilled if [</>J = 0, which, coupled with Eq. 20, yields 
K m[0V((/>£((/>))I n = 0. Hence the condition necessary for the validity of Eq. 8 on 17 is compatible 
with the standard interface conditions of continuum mechanics. Owing to this argument, we 
ultimately impose 

K m [<£V(<££(0))] • n = on S(t), t G (0, T max ], 

because it is this interface condition which is really needed in our problem. 

In the sequel we will assume that the time evolution of the interface S(t), and consequently 
those of the sub-domains 17t(£), &>H(t), is given. In other words, in the subsequent a priori 
estimates we will disregard the explicit coupling of their dynamics with the ones of the cells. We 
defer to a forthcoming work a more detailed investigation of such additional issues of the problem. 
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3.3 Initial conditions 



The initial distributions of cells and nutrient in il are described by two functions <j>o, Co : ft — > R, 
which should not exceed the expected physiological ranges. Therefore: 

(H8) < 0o < max , < c < c b in £1 

In particular, as observed in Sect. 3.1, nutrient concentration is expected to stay below the average 
physiological value c& because, in the present context, blood is the only source of nutrient, which 
is then consumed by cells. 



4 Notations and theoretic background 

In this section we prepare to tackle the analysis of the above-discussed models. We fix the main 
notations and quickly recall some essential technical material. 

Domain. The spatial domain is an open and bounded set !lcl <! (physically d = 1, 2, 3) with 
smooth boundary d£t. We use a for the coordinate along Oft and da for the (d — l)-dimcnsional 
Hausdorff measure in R d . The time interval of interest is [0, T max ], with finite final time T max > 0. 
We denote by Q Tmax the cylinder O x (0, T max ] C R d+1 . 



Functions. In general, we regard (j>, c as functions parameterized by time. Particularly, when 
we want to emphasize the dependence on t we write (j)(t), c(t) for the functions of x defined as 
Mt))(*) = *), (c(t))(x) = c(t, x). 

We will occasionally denote the time derivative by the subscript t for short (e.g., <j> t — d t (f> 7 
c t = d t c). 

We write Lip(u) for the Lipschitz constant of a function u, and u + , u~ for its positive and 
negative part: u + := max{u, 0}, u~ :— max{0, — u}. 

The indicator function of a set A is 1a, i-e., 1a (x) = 1 if x £ A, 1a{x) = if x £ A. 

Function spaces. L 2 (Q) is the usual Hilbert space of square-integrable functions in fi, endowed 
with the norm 

1/2 

(x)\ 2 dx 

We will also use a weaker L 2 norm denoted by ||| • |||o (cf. Appendix A). 

i? 1 (fi) is the Sobolev space of L 2 functions with square-integrable (weak) derivatives in il, 
endowed with the norm 

IH|i:= (IMI^ + llV-lloO 172 . 
For u e _ff 1 (fi), Stampacchia's Theorem asserts that u + , u~ G _ff 1 (fi) as well, with Vu + = 
Vul {u>0 }, Vu~ = -V?il {u<0} . 

Hq B (il) is the space of H 1 functions whose trace vanishes on B C d£l. The L 2 norm of the 
trace along B is 

1/2 



|o,B := {j h| 2 ^j 



We will deal, in particular, with the case B = c^Q. 

L 2 (0, T max ; H 1 (fi)) is the space of functions of t, valued in i7 1 (fi), which are square-integrable 
in the interval [0, T max ], endowed with the norm 
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The spaces L 2 (0, T max ; L 2 (Q)), L 2 (0, T max ; Hq g n (Q)) are defined analogously, and their respec- 
tive norms denoted similarly. In particular, in the former we will use the norm 




I • HI o being defined in Appendix A. 

We introduce the following shorthand notations: 

. V Tmax := L 2 (0, T max ; H 1 ^)) x L 2 (0, T max ; H^Sl)) 

• V := H 1 ^) x H 1 ^) 

• V f := ffoUn^) 

• for the dual space of Vf . 

We use the symbol (•, •) for the duality pairing between Vf and Vf. Given u e L 2 (0, T max ; Vf) 
with u t e L 2 (0, T max ; V^), it results 

(ut(t), u(t)) = ~Ht)\\l 

We use the abbreviation "a.e." for properties which hold "almost everywhere" with respect to 
the Lebesgue measure. 

If / C R is an interval, C°(J) is the space of continuous functions in I, endowed with the norm 

II ' Hoc :=sup|- (x)\. 

xEl 



Estimates and constants. We write 

a < b to mean 3 C > : a < Cb, 

the constant C being understood to be independent of both a and b, when the specific value of C 
is unimportant. In this case, C may vary from line to line within the same computation without 
explicit notice. 



Inequalities. 

• Cauchy's: for all a, b e R, ab < + |^ with arbitrary e > 0. 

• Poincare's: if u € Hq B (Q) then ||u||o < ||Vu||o- The constant involved in this estimate, 
denoted by Cp, depends in general on Q and B. 

• Cauchy-Schwartz's: | J n u(x)v(x) dx\ < |ju||o|M|o for all u, v <G L 2 (il). 



5 The time-dependent problem 

In this section we consider the initial/boundary- value problem 



— - K m A$(0) 



dc 
di 



K m V$0) • n 



— DAc = 



0, 



Km [V$(</»)]-n = 
0(0) = 00, 



r(t, x, 0, c) 

Q(t, x, <j), c) 

-DVcn 

c 



T}(C - Cb) 
Cb 



c(0) - c 



m 



n r (t), n H (t) 



t€(0, T max ] 
in QT max 

on d b n x (0, T max ] 
on dffl x (0, T max ] 
on S(t), t G (0, T max ] 
in 51 



(21) 
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along with the series of hypotheses (H1)-(H8), and we look for estimates of nonnegativity, bound- 
edness, uniqueness, and continuous dependence on the data of the functions (f>, c. We assume 
that solutions exist to this problem, in the suitable sense specified below. Notice that, in 21, the 
interface condition has been conveniently rewritten in terms of V3>(^>) using Eq. 10. 

Definition 2 (Weak solutions for the time-dependent problem). A weak solution to Problem 21 
is a pair ((f), c) £ Vr max such that: 

(i) fa, a £ L 2 (0, T max ; V}) 

(ii) $(0) e L 2 (0, T max ; H\n)) 

(iii) <j> — fa, , c = c b on df£l x (0, T max ] in the trace sense 

(iv) (f)(0) = (f> £ L 2 (U), c(0) = c G L 2 (^) 
which satisfies 



(fa, vi) + (c t , v 2 ) 

+ / (/c m V$(0) ■ Vv 1 +DVc- Vv 2 ) dx + n / (c-c b )v 2 da 




for all v\, V2 £ V/ and a.e. i € [0, T max ]. 

5.1 Nonnegativity and boundedness of the solution 
Theorem 1. Any weak solution ((f), c) £ Vr max to Problem 21 satisfies 

< (f>(t, x) < </> max , < c(t, x) < c b for a.e. (x, t) £ Q Tt 



Proof. First we establish that (f>, c are a.e. nonncgative by showing <j> , c — a.e. in Q Tin 
Choosing Vi = (f)~(t), v 2 = c~(t) as test functions in Eq. 22 reveals 

-~{H-(t)\\l + \\c-(t)\\l)-n m [ &(-(f)-Wd?-\ 2 dx-D\\Vc-(t)\\l 



(t)\\l idb n - VCb J c da 



= E 

a=T, i 



J I E 7^^(-^)3a(c)^+<5(^) 2 -A Q ^(0)<z Q (- C -) C - \dx. 



Because of hypotheses (H2)-(H7), the right-hand side is nonnegative for a.e. t £ [0, T max ]. 
Integrating from to t < T max and using fa^ , = (hypothesis (H8)) we get then 



t 

H~(t)\\ 2 + \\c-(t)\\ 2 + 2 Km J J &(-cb-)\V4-\ 2 dxdT 

Q 

t t 

+ 2i?y \\Vc-(T)\\ldT + 2f 1 J \\c-(T)\\l dbn dT + 27 1 c b J c-da<{) 



o o d b n 
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for all t E [0, T max ], whence , c = a.e. in <2T max due to the nonnegativity of each term at 
the left-hand side (use hypothesis (HI) for the term containing $'). 

Next we claim (0 — max ) + = (c— q,) + = a.c. in C>T max , which amounts to < max , c < c&. 
Let := (0 — 0max) + , c :— (c — Cf,)+ for brevity. Taking vi — 0(f), v 2 = c(t) as test functions in 
Eq. 22 yields 



(lW*)llo + l|c(t)||o) +«m / $'(0 max + 0)|V0| 2 ^ + C||V5(t)||§ 



l_d 

Ydt 

+ V\\c(t)\\ld b n 

= E / ( E Ta/a^max + ^aW-^xnax + 0)0 

- X a h a ((/))q a (c b + c)dj dx, 

the right-hand side being this time nonpositive for a.e. t G [0, T max ]. Integrating in time and 
using now 0(0) = c(0) =0we obtain 

t 

H0(*)llo + l|c(t)||o + 2K m / / $'(0*+0)|V0| 2 dxdr 



o n 



2D J \\Vc(r)\\ 2 Q dr + 2r ] | ||c(r)||§ An dT<0 



for all t e [0, T max ], whence the claim follows by arguing like in the previous point. □ 

5.2 Uniqueness and continuous dependence on the initial data 

Theorem 2. Let (0j, q) G VT max , i = 1, 2, 6e £wo weafc solutions of Problem 21 corresponding to 
the initial conditions (0,,o, c, ; o) € £ 2 (0) x L 2 (0). TTien 

T max 

|||02-0l|||| ?i 2+ jf" | ($(0 2 )-$(0l))(02-0l)^di 



Q 



+ 



1 max 

l c 2- c illi? H i+ y Il(p2-ci)(t)||g iafc n* 



< 1102,0-01,0110 + 1^2,0-01,011^. (23) 

In particular, the solution corresponding to a given initial condition is unique. 

Proof. It is sufficient to prove the estimate 23, for then uniqueness easily follows out of it with 

0i,o = 02,o and c lfi = c 2 ,o- 
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For all vi, v 2 <E Vf, the difference (fa — fa, c 2 — c\) of the two given solutions satisfies 

((02 - <f>l)t, Vl) + ((C 2 - Ci) t , W 2 ) 

+ y (K m V($(0 2 ) - $(0!)) • V«i + DV(c 2 - ci) • Vv 2 ) dx 



+ V J(c2-ci)v 2 da 



E /( E 7a{[/a(^2) - miWate) + /a (01 MM ~ <(d)]}«l 

_T U •> V d 

<S(0 2 - 0l)«l - A a (/l a (0 2 ) ~ /la(0l))9a(C2)«2 



a=T > H n a (t) W P< 



+ A Q ft, Q (0i)(g a (c 2 ) - g Q (ci))w 2 dx, (24) 



whence, choosing the test functions «i = ^(0 2 — fa), v 2 = c 2 — c\ and using Eq. 42 (cf. Ap- 
pendix A), we rewrite the left-hand side as 



l.h.s of (24) = 1| (1(02 - + life - Cl )(t)||g) 



+ 



I {$>(fa)-${fa)){fa-fa)dx 



+ ^||V(c 2 -c 1 )(t)||2+r ? ||( C2 -c 1 )(i)||^ n . (25) 

As for the right-hand side, Cauchy-Schwartz's inequality for the standard inner product in 
L 2 (f2), along with the boundedncss of /£, . . . , q a in the ranges of 0, c (recall Theorem 1), allows 
us to bound it from above as 

r.h.s of (24) < ]T 7S(llflalloo||/^)-/S(0i)||o 

a=T,H 
iy—p. d 



+ ll/alloo||^(C2) - <£(C1)|| )|W2 - 0l)(t)|| O 

+ \\h a {fa) - /la(0l)||ol|(c 2 - Ci)(t)||o, 



a=T, H 



where we have further used that — A a /i a (0i)(g a (c 2 ) — «7 a ( c i))( c 2 — Ci) < a.e. in QT max because 
(7 Q is nondecreasing (hypothesis (H7.2)). Since f£, h a are $-Lipschitz continuous, it results 

\\f v a {<h) - f v M\\l \\K{fa) - K{fa)\\l < f (Hfa) Hfa)) (h fa) dx, 

Combing this with the Lipschitz continuity of the g^s (hypothesis (H4)), Cauchy's inequality, and 
the fact that \\^ ■ ||o ^ ||| ■ |||o (cf. Appendix A), after some algebraic manipulations we arrive at 

r.h.s. of (24) < e J ($(0 2 ) - $(00) (0 2 - fa) dx 

+ Q + 7) (K& - + life - c.mWl) (26) 
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where e > is arbitrary. From Eqs. 25, 26 we deduce then that there exists C > such that 

~ (IllOfe -<Pl)(m\ 2 + \\(C2-C 1 )(t)\\l) 



+ (« m - eC) J (*(&) - $(0!)) (0 2 - fa) + D||V(c 2 - Cl )(t)\\ 2 

+ V\\(C2 dXOIIiUn < C Q + 1) (1(^2 - 0i)(t)|g + IKca - d)(t)||g) 

for a.e. i S [0, T max ]. Particularly, it is possible to choose e so small that K m — eC > 0. Multiplying 
both sides by e- 2C '\ C := C(l/2 + 1/e), yields 

~ (^ 2C ' 4 (1(^2 - + 11(02 - Cl)(t)||g)) 

+ e- 2C '*( Km -eC) y (*(02)-$(0i))(02-0i)da: 



n 

+ e-^'t (D\\V(c 2 Cl )(t)\\l + 7,11(02 - d)(t)||g An ) < 
whence, integrating from to t < T max and considering that e 2C * < e 2C Tmax for all < t < r max , 

^(l(^-0l)(*)io + ll(C2-Cl)(t)||g) 

t 

+ (n m -eC) J J ($(fa)-$(fa))(fa-fa)dxdr 
o n 

t t 
+ 2D J ||V(c 2 -ci)(r)||gdr + 2T ? f ||(c 2 - ci)(r)||g An rfr 


< l g 2C'T max (|||^ 2 o _ ^ i o |2 + || C2 o _ Ci o ||2) 

for all i e [0, Tmax]- At this point it suffices to observe that each term at the left-hand side, being 
nonnegative, is singularly bounded from above by the right-hand side. Integrating the first twos 
on [0, T max ] and evaluating the remaining ones for t = T max gives the thesis. □ 



6 The stationary problem 

In this section we turn our attention to the stationary problem 

-K m A$(fa = T(x,ct>,c) 
-DAc = Q(x, cj), c) 

* K m V<&(0) • n = 0, -DVc • n = r?(c - c b ) 

<j> = fa, c = c b 

{ K ro |[V$(0)] ■ n - 

which describes the equilibrium configurations of the model for large times. The asymptotic 
stability of constant steady states for a reaction-diffusion system sharing some analogies with 
Problem 27 has been addressed in [23]. 

Hcuristically, the solution of Problem 27 is what the solution of the time-dependent Problem 21 
tends to for t — > +oo. Making this limit rigorous with the appropriate concept of convergence is 
beyond the scope of this work, therefore we will be satisfied with the above intuitive interpretation. 

In Problem 27 the state variables depend on space only: <j> — fax), c = c(x). The sub-domains 
Qt, &h are fixed, their interface being S. 

We assume that solutions exist to Problem 27 in the following sense: 



in fir, flu 
in Q 

on 9b O 
on <9/f2 
on S 



(27) 
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Definition 3 (Weak solutions to the stationary problem). A weak solution to Problem 27 is a 
pair (</>, c) G V such that: 

(i) <&(<£) Gtf 1 ^) 

(ii) = , c = Cb on in the trace sense 
which satisfies 

J (K m V$(0) • Vui + DVc • Vu 2 ) dx + n J (c-c h )v 2 da 

= E / I E 7a/a(«(c)«l-^l-Aafea(0)?a(c)«2 I dx (28) 

for all v\, V2 G Vf. 

6.1 Nonnegativity and boundedness of the solution 
Theorem 3. Any weak solution (0, c) G V to Problem 27 satisfies 

< 4>{x) < 0max, < c(x) < q, for a. e. x G 17. 

Proof. In order to prove that 0, c > a.e. in 17 we choose v\ — <f>~ , v 2 = c~ in Eq. 28 and, 
similarly to the proof of Theorem 1, we compute: 

-J {n m ^'{-^-)\V4>-\ 2 + D\Vc-\ 2 ) dx + n J (c-c^c" da 
n d b n 

= E /( E ^/ a "(-fK(c)f + 5(^-) 2 - \ a K{4>)q a {c-)c- I dx. 

The right-hand side being nonnegative in view of hypotheses (H2)-(H7), it results 

-K m y$'(-</' _ )|V^| 2 dx-L>||Vc-||2-r ? ||c-||2 9b0 -77c b J c~ da>0 7 
n d b n 

whence, considering that each term at the left-hand side is nonpositive, we deduce immediately 
cT = a.e. in 17, and obtain furthermore <&'(— <f>~)\ V</>~ | 2 = a.e. in 17. This means either 
<£>'(— 4>~) — 0, which implies <p~ — a.e. in 17 because $' vanishes at most in zero, or |V</>~| 2 = 0, 
which yields <j>~ = a.e. in 17 as well due to <p~ G V) (use Poincare's inequality). 

Next we prove that </> < ma x, c < q, a.e. in 17. Set <j> := (0 — </> max ) + , c := (c — q,)+ and 
choose wi = </>, «2 = c as test functions in Eq. 28 to find 

n m J $'(<?W + 0)|V0| 2 dx + D\\Vc\\l + r/||S|| 2 9b0 
n 

= E /( E 7£/aWm« + #7a(c)0 

-6(<f> 

mux 

+ 0)0 - \ a h a (4>)q a (c b + c)c dx. 
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Arguing like in Theorem 1, we conclude that the right-hand side of the above equation is nonpos- 
itive, whereas each term at the left-hand side is nonnegative. Therefore we have c = and also 
&'{<f>max + 0)|V0| 2 = a.c. in ft. Since $'(0 ma x + 0) > 0, this implies |V0| 2 = a.e. in 0, i.e., 
= a.e. in as well thanks to <G Vf. □ 

6.2 Uniqueness of the solution 

Theorem 4. There exists a constant C > 0, depending only on the coefficients X a and on the 
functions f£, g v a , h a , q a , such that if C is sufficiently small then Problem 27 admits at most one 
weak solution (0, c) e V. 

Proof. Let (0j, a) £ V, i = 1, 2, be two solutions, then for all v\, V2 <E Vf the difference ((f>2 — 
(f>i, C2 — ci) satishes 



J (« m V($(0 2 ) - $Oi)) • V«i + DV(c 2 - ci) • Vv 2 ) 
n 

+ ^/ y (C2 - ci)«2 da 

= E /( E 7^{[/^(02)-/^l)]^(C2) + ^(0l)^(C2)-^(ci)]}t;i 
a=T, ifj^ \"=P, d 

- (5(02 - 0i)«i - K{K{4>2) - h a ((t>i))q a (c 2 ) 
+ A a /i a (0i)(g a (c 2 ) - q a {ci))v 2 dx. 



We choose v\ — ^(0 2 — 0i) (cf. Appendix A), u 2 = c 2 — c\, and mimic the computations of the 
proof of Theorem 2 (using in particular Cauchy's inequality at the right-hand side with e = 1) to 
deduce 



Km /(*(&) - *(0l))(02 - 0l) dx + D\\V(C 2 - Ci) 



n 

Hence there exists C > such that 



+ »?l|C2-C 1 || 2 An + ^|||02-0l|||o 

< I (^-^iVifo-^dx + ifo-Ml + llci-c^l. 



(n m -C) J ($(0 2 )-$(0i))(0 2 -0i)cte + (<5-C)|||0 



2 - 01 Wo 



52 



D-CC% 
p 



where we further applied Poincare's inequality to c 2 — c\ G Vf. Uniqueness of the solution follows 
from this relationship as long as the coefficients of all terms are positive, which entails 

C < min j/c m , 6, . □ 
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Remark 1. For the sake of definitcness, we record that a possible constant C for Theorem 4 is 

E 7aLip$COI.<Ooo + A Q Li P$ (MMI 




E E 7^ll^llooLip 2 (.^) + A Q || (ZQ || oc , 

a=T,H \i/=p, d 

c% E 7^ (Halloo + 

a—T. H ) 



f=p, d 

6.3 Existence of the solution 

We complete our analysis of model 27 by outlining the theory of the existence of solutions. We 
confine ourselves to the one-dimensional setting, taking as reference domain the dimensionless 
interval / = (0, 1). In particular, x — will be the vascular boundary and x = 1 the far boundary. 

The case d = 1 allows us to rely on two basic tools, which are not available in higher dimensions: 
on the one hand the Sobolev embedding C°(J) C H 1 ^), on the other hand Morrey's inequality 
||w||oo < |M!i for u e Extending the theory to the case d > 1 is likely to require partly 

different tools, which is at present beyond the scope of the work. 

The one-dimensional problem is written as 

-«m*(0)xs - r(z, 0, c) in (0, S), (S, 1) 

-Dc xx = Q(x, 4>, c) in / 



K m $(cj>) x (0) = 0, -Dc x (0) = V {c(0)-c b ) 

<Hi) - <t>*, c(i) - c b 



(29) 



where S € / is the location of the point interface between tumor and host cells. In particular, 
CIt = (0, S) and = (S, 1). By adapting Definition 3 to the present context, a weak solution 
to Problem 29 is a pair ((/>, c) e V, such that € H 1 ^), 0(1) = 0*, c(l) = c b , which satisfies 



(K m $(0) x t>ix + Dc x v 2x ) dx + rj (c(0) - c b ) v 2 (0) 



o 



= E / E 7:/ a "(^(c)«i-%i-A a ft tt (^)« a (c) 1)2 & (30) 

for all «i, t)2 € V/. 

For the subsequent theory, it is useful to introduce the inverse <I> _1 of the constitutive function. 
Owing to hypothesis (HI), $ _1 is continuous and strictly increasing, with <& _1 (0) = 0, <& _1 (s) < 
for s < 0, and <I> _1 (s) > for s > 0. It is also smooth on (— oo, 0) and (0, +oo), however 
linis^o^ -1 )'^) = +°o because of the degeneracy of at the origin. 

In order to prove the existence of stationary solutions we resort to a splitting method, which 
consists in approaching the two differential equations of Problem 29 separately, assuming that the 
main unknown is either <j> or c and that the other function is known. 

Theorem 5. Assume h a (0) = 0. There exists a constant C > 0, depending only on the parameters 
K m, 7a; ^max and on the functions $ ; fa, g v a , such that ifC is sufficiently small then Problem 29 
admits a weak solution (<ft, c) G V. 
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Proof. We preliminarily define the sets 

V := {/ G L 2 {I) : < / < max a.e. in 1} 
U := {/ e i 2 (I) : < / < c 6 a.e. in J}. 

Let us begin by considering the problem 

-Dc xx = Q(x, ip, c) in / 

Dc x (0) = n(c(0)-Cb) (31) 
c(L) = c fc , 

where (5 e V is given. We associate with it an auxiliary problem in which the function q a is 
replaced by q a — q a l[o, +<»)• The corresponding weak formulation is obtained from Eq. 30 by 
letting vi = and writing g Q in place of g a : find c G H 1 {I) 1 with c(l) = q>, such that 

l 

D c x v x dx + r?(c(0) - Cb)v(Q) = - ^ / \ a h a (<p)q a (c)v dx (32) 



a=T, ff, 



for all f G V/. By introducing the antiderivative Q a of g Q vanishing in zero, we can view Eq. 32 
as the Euler-Lagrange equation for the functional 



l 



Jl(c) = y J C 2 x dx+^ ( C (0) - C 6 ) 2 + J2 j X Mf)Qa(c) dx 

over the class of admissible functions A\ = {c G H 1 ^) : c(l) = c^}. Thus we can seek our solution 
c as a minimizing point of J\ on A\ . 

Since q Q (s) = for s < and q a (s) = q a (s) > for s > 0, we have Q Q (s) > for all set. 
Using further that A Q and h a are nonnegative we obtain Ji(c) > D/2||c x ||o, which implies that 
Ji is coercive. Therefore any minimizing sequence {ck}^ =1 Q A\ is bounded in i? 1 (I) and, upon 
passing to a subsequence, we can assume that it converges weakly to some c G ff 1 (0, L). But 
Cfe — Cb G V/ all fc and Vf is a weakly closed subspace of i? 1 (/) (in view of Mazur's Theorem, as 
it is closed), thus we deduce more precisely c— q, G Vf, i.e., c(l) = q, and ultimately c G Ai. 

Considering that J\ is of the form J Q L\{c x , c, x) dx for the Lagrangian 

Li(p, z, x) = yp 2 - r](z - c b )p+ x »(x)h a ((p(x))Q a (z)l na (x), 

a=T, H 

which is smooth and convex in p for all z G R and all x G /, we deduce that J\ is sequentially 
weakly lower semicontinuous on i? 1 (7). Thus c is a minimizing point of Ji, i.e., a solution to our 
auxiliary problem. 

Mimicking the computations of the proof of Theorem 3 with v\ = reveals that, for any fixed 
<p G V, all solutions to the auxiliary problem range in [0, c&]. Hence we conclude < c(x) < Cb 
for all x € I, and consequently that c solves also Problem 31 because the latter and the auxiliary 
problem coincide for c G [0, c b \- Notice that c G U. 

We show now that the solutions to Problem 31 depend continuously on ip. Let ci, C2 be two 
solutions corresponding to ipi, <p2 € V, respectively, then for all v G Vf the difference C2 — ci 
satisfies 

l 

D /(ca - ci) x « a + 77 (c 2 (0) - ci(0)) v(0) 

= - 51 / A «{(/la(<^2) - MVl))9a(c 2 ) + /la(yi)(9a(c 2 ) ~ q a {ci))}vdx. 

T 1 IT *^ 
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We chooser = c 2 — c\ and observe that — \ a h a ((pi)(q a (c 2 ) — q a (ci))(c 2 — Ci) < in fl a (hypothesis 
(H7.2)), whence 

£>||(C2 - ci^Hg + (ca(0) - d(0)) 2 



< 



E / ^aq a {c2)\h a {^2) - K{ipi)\ ■ |C 2 - Cl\d,X 

\ E A «ll9«l|oo ( -\\h a {<fi2) - M^i)llo + e\\c 2 - d\\l J 



< 

2 



Now we recall, from hypothesis (H6), that h a is $-Lipschitz continuous in [0, </>*], whence 

l 

\\h a (<p 2 ) - M<£i)llo < Lip*(/i«) y (*(V2) - *(</?i)) (</>2 - </?i) cfe 

o 

^ ll<P2 - ^lllo- 

Thus the previous computation can be continued by asserting that there exists C > such that 
D — eCCp . 2 2 C,| 1,2 



1 + C 



2 



C2-c 1 ||f + r ? (c 2 (0)-c 1 (0)r < -H^-ViHo. (33) 
y p e 

where we also applied Poincare's inequality to c 2 — c x e Vf. Choosing e > so small that 
D — eCCp > 0, we get from Eq. 33 the desired continuity estimate. In particular, for ipi = ip 2 we 
obtain the uniqueness of the solution to Problem 31. 

The foregoing results enable us to define the operator :V^U such that S^i (ip) = c. From 
the continuity estimate 33 we deduce that Sf\ is Lipschitz continuous on V and, with a little more 
work, that it is also compact. To see this, we observe first of all that the assumption h a (Q) = 
implies ^i(O) — ci, (i.e., the unique solution to Problem 31 for ip = is c — q,), then we choose 
ipi = in Eq. 33 and drop the subindex 2 to obtain 

ll^i(¥>)-<%||?<|M|g. 

We take now {<Pk}kLi C V and notice that, in view of the latter estimate, the sequence {,5^\{ipk) — 
Q>}fcLi is bounded in iJ 1 (7). Owing to Rellich's Theorem, we can therefore assume, upon pass- 
ing to a subsequence, that .y\{<pk) — Cb converges in L 2 (I) as k — > oo, i.e., that the sequence 
{<yi{<Pk)}kLi is convergent, which proves the compactness of S?\. 



We turn now our attention to the problem 

x, 9) in (0, S), (S, 1) 

(34) 



-Km${<!>) xx = T(x, 0, 9) in (0, S), (5, 1) 

Km*Wx(0) = 

0(1) = </>» 

K m M4>) X i = 0, 



where e U is given. Again, we associate with it an auxiliary problem in which the functions 
fP, f* are replaced by fg = /£l[o,<£ max ], /„ = /q1[o,+k,), respectively. The weak formulation 
is recovered from Eq. 30 by letting v 2 — and substituting conveniently the functions at the 
right-hand side: find <p € H 1 ^), with $(0) e H 1 ^) and 0(1) = 0*, such that 



, / *(4>) xVx dx = E / ( E <r a WM - sA vdx 

J T 1 U J \ ,1 I 



1 

for all v € V/. We set u := $(0), whence = $ _1 (w), so that this equation becomes 

n m Ju x v x dx=Y J J[ E ^(^r'j^^-^wLx (35) 
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for all v £ Vf. If we introduce the antiderivatives F£, ^ of o 4> x , $ 1 vanishing in zero, we 
can regard Eq. 35 as the Euler-Lagrange equation for the functional 



1 



■h(u) = ^ fuldx- ]T ( £ -f a FZ(u)g v a (0) <**(«) d* 
o Q=T . H n a \"=J>. d / 

over the class of admissible functions _4 2 = { u G : u(l) = <f>(0*)}. Thus, again we can look 

for our solution u as a minimizing point of J 2 on A 2 • 

Notice that Pp(s) < FP($(0 max )), and that F*(s), *(s) > for all s€l. Therefore, recalling 
further that 7„ < (hypothesis (H2)), we have 

Mu)>^\\u x \\l- ]T i p Jg p a \\ooF^(cb max m a \. 

a=T, H 

J 2 is thus coercive, hence any minimizing sequence {uk}^ =1 C _4 2 converges weakly (up to possibly 
passing to subsequences) to some u £ H 1 ^). Since Uk — $(</>*) £ Vf and Vf is weakly closed 
in -H' 1 (/), it results u — $(0*) £ Vf, that is u £ Ai- In addition, J 2 is in turn of the form 
L 2 (u x , u, x) dx for the Lagrangian 

L 2 (p, z, x) = ^p 2 - ^ 7^(*)<?»)) - ln a (s), 

a=T,H \f=p, d J 

which is smooth and convex in p for each z £ R, x £ /, hence J 2 is sequentially weakly lower 
semicontinuous on i? 1 (/). It follows that u is a minimizing point for J 2 on .4 2 , and consequently 
<f> := is a weak solution to our auxiliary problem. 

Mimic now the computations of the proof of Theorem 3 with v 2 = to obtain that, for any 
hxed 9 £lA, all solutions to the auxiliary problem range in [0, max ], whence < (f> < </> max in I. 
But the auxiliary problem and Problem 34 coincide for 4> £ [0, max ], hence ultimately we have 
found a solution £ V to Problem 34. 

Next we show that, by introducing suitable constraints on the parameters, we can guarantee 
that </> be strictly positive in I. Let us pick v\ — $(</>) — $(0*), v 2 = as test functions in Eq. 30 
to discover 

l 

K m f $(0), ($(0) - H4>*)) x dx= J2 I 7S/a(«W W) - dx 

u—p, d 
1 



-6 y </>($(</>) -$(</>„)) dx. 



Noting that §{<p) x — ($(</>) — ®(4>*)) x at the left-hand side and using the boundedness of /£, 
at the right-hand side, we estimate 

||(*W-*(0.))*llo^ll*W-*(^)ll=o- 

In addition, owing to Poincare's and Morrey's inequalities, 

ik$(0) - $(^)),n§ > \\m > wm $(- ■ < 



I oo ' 



hence finally there exists C\ > such that || $(</>) — $(</>*) ||oo < Ci, which indicates that $(0(x)) > 
$(0*) — C\ for all x £ I. For dehnitcncss, we report the explicit expression of a possible constant 
Ci: ' 

/ \ 



c 1 = l±^ 



E 7SII/alloo|^lloo|fia, 



a=T, H 
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We fix now e £ (0, </>*) and observe that <f> > e if and only if <£>(</>) > 3>(e), thus we can guarantee 
that 4> be strictly positive in I if we require $(</>*) — C\ > $(e), which implies the constraint 

Ci < - $(e). (36) 

Given this, any solution <f) £ V to Problem 34 satisfies < e < 4> < '/'max in /. 

Finally we assert that, under condition 36, solutions to Problem 34 depend continuously on 
6 £ U in the norm || • 1 1 o - For this, let <f>i, <j>2 be two solutions corresponding to 6\, 82 £ U, 
respectively, then for all v £ V/ their difference <f)2 — (f>i solves 

1 1 

Km J(^((f>2) ~ $(4>l))xV x dx + S j {(f>2 - <f>l)vdx 


= E /7«{(/«(02)-/:(^i))^(e 2 ) + /:(^ 1 )(^(0 2 )-^(0 1 ))}«^. 

v—p 7 d 

We choose v = &(4>2 — 4>i) (cf- Appendix A) and, mimicking the computations of Theorem 4, we 
find that there exists C2 > such that 

1 

(Km-C 2 ) I ($(02)-*(0l))(02-0l)d!C + (*-C2)||^-0l||g<||fl2-(9l||§ 


(for the sake of completeness, we point out that the constant C 2 is the same as the one appearing 
in Theorem 4, cf. also Remark 1). Assume 

C 2 < min{« m , <5}, 

then, since |$(</> 2 ) - $(0i)| > (min se[ej 0max] $'(s)) l<fc - </>i|, it follows 

min $'( S )(« m - C 2 )||0 2 -faWl + tf- C 2 )U2 - 0i|g < ||(9 2 - 041 (37) 

which yields the desired continuity estimate, together with uniqueness of the solution to Problem 34 
when 6\ = 9 2 - 

Define now C := max{Ci, C 2 } and impose 

C < min{$(0») - $(e), K m , 5}, 

then Problem 34 admits a unique solution <fi £ V for any given 9 £U. Consequently, we are in a 
position to define the operator 5^2 : U — > V such that ^(6*) = 0, which is Lipschitz continuous 
on W in view of Eq. 37. 

At last, we come back to the full Problem 29 in this way: we construct by composition the 
operator := =5*2 o 5?\ : V — >• V such that =5^(<p) = </>. Since «5*i is continuous and compact and 
■5*2 is continuous, ^ is in turn continuous and compact; moreover, V is convex and closed in L 2 {I) 
(to see the latter property, use that convergence in L 2 (I) implies pointwise convergence a.e. in / 
upon passing to subsequences). Schauder's Fixed Point Theorem implies then that 5? has a fixed 
point (f> £ V, hence the pair (<f), c — (</>)) is a weak solution to Problem 29 and we are done. □ 

Remark 2. If, in addition to the hypotheses of Theorem 5, also the hypotheses of Theorem 4 hold 
true then the solution to Problem 29 is unique. 
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7 Possible developments 



In this paper we have addressed the mathematical formulation of initial and boundary-value 
problems for multiphase models of tumor growth, deduced from the framework developed in [32] . 
We have performed a qualitative analysis of both the time-dependent and the time-independent 
problems, mainly by means of L 2 -H 1 a priori estimates, establishing nonnegativity, boundcdness, 
uniqueness, and continuous dependence of the solution on the initial data. In the one-dimensional 
time-independent case we have also obtained the existence of the solution. 

The analytical techniques used here may be profitably exploited to approach more advanced 
multiphase models, also fitting the framework presented in [32], which incorporate a more ac- 
curate description of the interactions between the cells and the extracellular matrix. Based on 
phcnomenological laboratory observations [5, 18, 33], they take into account the adhesion of the 
former to the latter by relating the cell-matrix stress m am (a component of the overall external 
stress m a included in Eq. 2) to the cell-matrix relative velocity v m — v Q as 

v m - v a = K am (l - Smn ) m Qm . (38) 



This formula says that if the magnitude of the stress m am is below some critical threshold £ am > 
then there is no relative motion between the cells and the matrix, that is the former remain 
anchored to the latter. Conversely, if |m Qm | is above the threshold S Qm then the interaction 
stress m am is proportional to the relative velocity v m — v a , thus recovering a more classical 
viscous friction which, in particular, means that cells slide on the matrix with their own velocity. 
If Eq. 38 is used, with the additional assumption of motionless matrix (v m = 0), then the equations 
ruling cell dynamics take the form 

^jr - V ■ [4> a s a {(f>T, 4> H , |V^|)K arn v(0S(^))] - r a , (39) 

where 

\ + 



translates the adhesion mechanisms discussed above. In particular, the velocity v Q of the cells is 

v« = -y a K arn v((f>^)), 

hence if |V(0£(</>))| < S Qm then v a = because = (recall that (f> a < <p by definition) and 
the cells stay attached to the matrix, while if |V(<^E(<£))| > S Qm the cells might detach from the 
matrix since one may have J? a > 0. 

Equation 39 can be regarded as a refined version of Eq. 3, which would be interesting to study 
in view of its physical significance, possibly adapting the techniques illustrated in this paper. 
Notice indeed that setting S Qm = for both a = T and a = H, which amounts to assuming a 
purely viscous friction between the cells and the matrix without attachment /detachment, reduces 
Eq. 39 to Eq. 3, hence the latter turns out to be a particular case of the former. Additional 
mathematical difficulties need however to be overcome, especially the harder degeneracy of the 
differential operator in space caused by the term J? a . 
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A The Poisson solution operator 



In this appendix we introduce a basic tool, inspired by [24, 29], useful to handle the nonlincarity 
$ofEq. 11. 

Consider the linear elliptic problem 



— Au = f in ft 
Vm • n = on <9hft 
u = on <9/ft 



(40) 



and assume <9/ft ^ (see Appendix B for the case dffl = 0), / G L 2 (ft). The classical weak 
formulation is: 

find u € Vf such that 



J Vti -\7vdx = Jfv dx, V v G V) , 



(41) 



II o - 



n n 

which, owing to Lax-Milgram theory, yields a unique weak solution u with ||Vit||o < 

Let us introduce the linear operator & : L 2 (ft) — > L 2 (ft), termed the Poisson solution operator, 
which associates with / the solution u = &f to Problem 41 above. is bounded, symmetric, 
and positive definite. Therefore, besides the standard inner product (/, g) = J n f(x)g(x) dx, we 
can endow L 2 (ft) with an inner product induced by ((/, gj) :— (£Pf, g), with corresponding 
norm := ((/, /)) = (&>f, /). Using & we rewrite Eq. 41 as 



J V0>f -Vvdx = J fvdx, V/ e L 2 (Q), v e Vf, 



(42) 



whence, letting w = ^/ G V/, we see that |||/|||o = HV^/Ho- Thus the previous a priori estimate 
on u implies |||/|||o ^ ll/llo- On the other hand, in view of Poincare's inequality, we also have 
IWIIo < llV^/llo HII/lllo. 

Given u G L 2 (0, T max ; Vf) with u t G L 2 (0, T max ; Vf), the duality pairing between u t (t) G 
and ^u(t) G V/ yields, for a.e. t, 

( Ut (t), = ~K*)Bo- 



B Problems with no far boundary 

The aim of this appendix is to extend the theory developed in Sect. 5 to problems in which there 
is actually no far boundary within the host environment, that is the whole tissue is surrounded by 
blood vessels (<9ft = c%ft, <9/ft = 0). For instance, this situation arises when studying the confine- 
ment of a tumor mass within a healthy tissue well supplied with blood by a nearby vasculature. 
In this case we can imagine that ft-r(i) C ft at all times, with dfl = dfln(t) = <9fcft (Fig. 5). From 
the mathematical point of view, the main point is the disappearance of the Dirichlet boundary 
conditions, which forces one to modify the Poisson solution operator & introduced in Appendix A 
in order to deal with the new problem for <j>. 

Specifically, in place of Problem 40 we consider (cf. also [24]) 

' -Au = /-(/) in ft 
< Vu • n = on dft 

^ («) = (/) in ft 

for / G £ 2 (ft), where (•) denotes average on ft: 

(/> :=^jf(x)dx. 
n 
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Figure 5: The geometrical configuration of the domain ft in the problem of a tumor mass sur- 
rounded by host tissue with no far boundary. 



The weak formulation is: 

find u G H 1 ^!) such that 
Jvu-Vvdx + \n\(u)(v) = Jfvdx, Vveff 1 ^), ( 43 ) 



whence, owing to Lax-Milgram theory, we get a unique weak solution fulfilling ||m||i < 
Setting u — & f ', we observe from Eq. 43 that the operator satisfies now 

Jv&f-Vvdx + \Sl\{&f){v)= J fvdx, V/ e L 2 (0), v e H 1 ^), 

therefore the inner product ((•, •)) and the induced norm ||| • |||q take the following forms: 

((/, 9)) = J V5»/ • V&gdx + \Sl\(&f){&g) 
n 

Ill/Ill 2 , = \\v&f\\l + \n\(&ff^ ||^/||?. 

This essentially affects the estimates of Theorem 2 when dealing with 
J V($(0 2 ) - • V^(0 2 - ^i) di 

($(&) - $(0l))(02 - 0l) - |Q|<$(^ 2 ) - $(0l)><5»(02 - 0l)>, 



10- 



because of the second term at the right-hand side which must be estimated. First we use Cauchy's 
inequality to find 

- $(</>!)>( ^(02 - 0l)) < | - *(0l)) 2 + ^(^(02 - ^l)) 2 - 
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Next we employ Cauchy-Schwartz's inequality and Lipschitz continuity of $ to discover 

- $(<M) 2 < ±- J - $(^i)| 2 dx 
n 

< J m<h) - *(0i))(02 - ^, 
w 2 - 4>i){t)? < T^wnfc 0i)(*)iig < ci(0 2 - 0x)(t)i!ie, 



whence 



y v(*(&) - $(^i)) • - 00 ^ 



Choosing e > so small that eLip($) < 2|f2|, the terms at the right-hand side can be finally 
incorporated into the similar ones already present in the proof of Theorem 2. 
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